Phase diagram of Josephson junction between s and s± superconductors in dirty limit 
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The s± state in which the order parameter has different signs in different bands is a leading can- 
didate for the superconducting state in the iron-based superconductors. We investigate a Josephson 
junction between s and s± superconductors within microscopic theory. Frustration, caused by in- 
teraction of the s-wave gap parameter with the opposite-sign gaps of the s± superconductor, leads 
to nontrivial phase diagram. When the partial Josephson coupling energy between the s-wave 
superconductor and one of the s± bands dominates, s-wave gap parameter aligns with the order pa- 
rameter in this band. In this case the partial Josephson energies have different signs corresponding 
to signs of the gap parameters. In the case of strong frustration, corresponding to almost complete 
compensation of the total Josephson energy, a nontrivial time-reversal-symmetry breaking (TRSB) 
state realizes. In this state all gap parameters become essentially complex. As a consequence, this 
state provides realization for so-called (^-junction with finite phase difference in the ground state. 
The width of the TRSB state region is determined by the second harmonic in Josephson current, 
oc sin(2</!>), which appears in the second order with respect to the boundary transparency. Using 
the microscopic theory, we establish a range of parameters where different states are realized. Our 
analysis shows insufficiency of the simple phenomenological approach for treatment of this problem. 



I. INTRODUCTION 

The discovery of superconducting iron pnictides and 
chalcogenides is one of the most remarkable recent 
achievements in the condensed-matter physics. A rapid 
progress in characterization of these materials and devel- 
opment of theoretical understanding has been reflected in 
several reviews^. The key feature of these semimetal- 
lic materials is the multiple-band structure, the Fermi 
surface is composed of several electron and hole pockets 
located near different points of the Brillouin zone. 

Superconductivity in the iron-based materials is likely 
to be unconventional. There is a theoretical consensus 
that the electron-phonon interaction is not strong enough 
to explain high transition temperatures^ In several theo- 
retical papers it was suggested that superconductivity is 
mediated by spin fluctuations leading to an unusual su- 
perconducting state in which the order parameter has op- 
posite signs in the electron and hole bands (s± state) 4^— 
Experimental verification of this theoretical proposal be- 
came one of the major challenges in the field. Probing 
the relative sign of the order parameter in different bands 
is not trivial and the structure of superconducting state 
has not been unambiguously established yet by experi- 
ment, even though several properties consistent with the 
s± state have been revealed. An extensive critical review 
of experiments both in favor and against the realization 
of the s± state in iron-based superconductors has been 
done recently in Refs. 0. Shortly, the main experiments 
supporting the s± state include 

• Observation by the inelastic neutron scattering of 
the resonant magnetic mode below the supercon- 
ducting transition temperature^ Such a mode is 
expected for the superconductors with the sign- 
changing order parameter. This mode was observed 
in almost all compounds and its frequency scales 
approximately proportional to the transition tem- 



perature. 

Microscopic coexistence of antiferromagnetism and 
superconductivity demonstrated in some com- 
pounds within a narrow doping range, most clearly 
in Ba[Fei_ a; Co a ;]2As2 i 10 ' 11 For the case of the con- 
ventional s ++ state in which the order parameter 
has the same sign in all bands, the spin-density 
wave (SDW) has a strong pair-breaking effect on 
the bands connected with the SDW ordering wave 
vector. Such direct pair breaking is absent if the 
order parameter in such bands has opposite signs 
meaning that the SDW is much more compatible 
with the s± state than with s ++ 
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• The magnetic field dependence of the quasiparticle 
interference peaks studied by scanning tunneling 
spectroscopy in FeSe^Tei- 
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On the other hand, discovery of the iron selenide com- 
pounds without hole band and with rather high transition 
temperatures, up to 30K— , questioned universality of the 
s± state for all iron-based superconductors. Also, it oc- 
curs that the iron-based superconductors are quite stable 
with respect to disorder. As for the s± state the inter- 
band scattering is pair breaking, stability with respect to 
disorder is frequently used as an argument against this 
state. Therefore the structure of the order parameter in 
the iron-based superconductors is an unresolved issue. 

One of the ways to probe unconventional supercon- 
ductivity is to study Josephson junctions and proximity 
effects with conventional superconductors. In the case of 
contact between s-wave and s± superconductors, frustra- 
tion, caused by interaction of the s-wave gap parameter 
with the opposite-sign gaps of s± superconductor leads 
to several anomalous features which were recognized and 
studied in several theoretical papers.— ~— For example, 
proximity with s± superconductor induces corrections to 
the density of states of s-wave superconductor which, in 
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FIG. 1. Left part: A Josephson junction between s and s± su- 
perconductors with the gap parameters for a general complex 
state. Right part: Generic phase diagram of such Josephson 
junction. 



principle, allow to identify the signs of the order param- 
eter in different bands of s± superconductor^ Particu- 
larly interesting is a possibility of a time-reversal sym- 
metry breaking (TRSB) state^^i^ in the parameter 
range where the partial Josephson coupling energies be- 
tween s-superconductor and different s± bands almost 
exactly compensate each other. Existing experiments 
on Josephson junctions between iron-based and conven- 
tional superconductors have been reviewed in Ref. |25| . No 
anomalous features, however, have been reported so far. 

For experimental realization of the TRSB state it is 
important to establish range of parameters where such 
state can be expected. Up to now this state was stud- 
ied using mostly phenomenological models which are not 
rigorously justified. The purpose of this paper is to de- 
velop microscopic description of transition between the 
aligned and TRSB states. The paper is organized as fol- 
lows. In Sec. HI] we present general consideration of the 
transition between the aligned and TRSB state in the 
region where the partial Josephson energies almost com- 
pensate each other. It is known that the width of the 
TRSB region is determined by the second harmonic of 
the Josephson current, which in our situation appears in 
the second order with respect to the coupling between 
superconductors. In Sec. Mil we present the microscopic 
equations and boundary conditions describing the con- 
tact between s-wave and s± superconductors in dirty 
limit considered in this paper. In Sec. IIVI we consider 
corrections to the Green's functions and gap parameters 
induced by the interface. Computation details of these 
corrections are presented in Appendix [X] The proximity- 
induced corrections to the Green's function determine 
the second harmonic in the Josephson current, which is 
considered in Sec. [V] We present both general formulas 
for different contributions to the second harmonic and 
simple analytical results for the most relevant limiting 
cases. We also reveal the dominating contribution to the 
second harmonic. Using these results we analyze in Sec. 
IVII the width of the TRSB region and its shrinking with 



increasing temperature. Finally, in Sec. IVIll we consider 
proximity-induced corrections to the density of states of 
the s-wave superconductor within the TRSB region. 

II. GENERAL CONSIDERATION OF THE 
TRANSITION REGION IN THE 
WEAK-COUPLING CASE 

We consider the Josephson junction between s-wave 
and two-band s± superconductors, see Fig. [TJ (left). 
In the weak-coupling limit this system is characterized 
by the partial Josephson coupling energies between s- 
superconductor and s± bands, Ej a with a being the 
band index. Typically, the s-wave gap parameter aligns 
along the s± gap with which it has larger coupling en- 
ergy. In this aligned state the partial Josephson cou- 
pling energies are positive and negative for the aligned 
and anti-aligned bands correspondingly. Nontrivial be- 
havior is expected in the case of strong frustration when 
the total Josephson energy is close to zero. This hap- 
pens when absolute values of the Josephson energies for 
the opposite-sign bands are close, \Eji\ ~ l-E^I- Phe- 
nomenologically, the phase diagram can be described by 
the model of the frustrated Josephson junction consid- 
ered in several papers^ which provides correct qualita- 
tive description. However, in general, this model does not 
describe the system quantitatively, because, ignoring the 
fermionic degrees of freedom in the s± superconductor, 
it does not treat correctly its interband energy. For the 
weak-coupling regime, however, the transitional region 
between the two aligned states can be treated following 
the same reasoning as for the transition between and 
7T junctions, see, e.g., Refs. HH andl26l In the vicinity of 
transition the linear approximation for the coupling be- 
tween the superconductors becomes insufficient and the 
total Josephson energy can be represented as 

£{<f) = (E n ~ E J2 )(1 - cos0) + -^-(1 - cos 2$, (1) 
where <f> is the phase difference between A s and Ai and 

(2) 

the term Ej appears in the second order with respect 
to the boundary transparency. This corresponds to the 
Josephson current 

3 (4) = Uji - 3J2) sin 4 + sin 20 (2) 

with jj a — (2-Kc/^o)Ej a . The intermediate TRSB state 
exists only if the sign of the second-harmonic is negative 

(2) (2) (2) 

Ej ,jj < 0. In this case in the region \ jji—jj2\ < 2jj 
the ground-state phase difference is given by 

cos0o-(jJi-X/2)/(2|jf ) |). (3) 

It smoothly transforms between and tt when the dif- 
ference jji — jj2 changes from 2\jj ' \ to — 2\jy '\. There- 
fore the TRSB state also provides realizations of so-called 
0-junction2i in which a finite phase difference exists in 
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ground state leading to several anomalous properties. In 

(2) (2) 

the case Ej ,j\ > the transition between the two 
aligned states is of the first order and the TRSB state is 
not realized. 

The simplest phenomenological description is the frus- 
trated Josephson junction model in which the tilt of the 
relative phase between two gap parameters of s± super- 
conductors is described by the energy £ i2 cos(#i — 6 2 )- In 
this model the amplitude of the second harmonic is given 
by 



- -jjEj/(2£ 12 ) 



J. j 



with j j = (jj! + JJ2 )/2 and Ej = (E, n + E J2 )/2. Our 
further microscopic analysis shows that this result is only 
valid for a special situation of very weak coupling between 
the bands of the s± superconductor. We will compute 
the second harmonic in general case within microscopic 
approach. 



III. EQUATIONS AND BOUNDARY 
CONDITIONS 



In this section we write down equations and boundary 
conditions for the simple microscopic model describing 
a "sandwich" , consisting of slabs of two-band s± super- 
conductor with thickness d± and a single-band s-wave 
superconductor with thickness d s , as shown on Fig. [1] 
We denote the bulk critical temperatures of the s-wave 
and s± superconductors as T c s and T c , respectively. The 
x = plane coincides with the interface between the su- 
perconductors. The main assumption of our description 
is that both superconductors are in dirty limit but the 
interband scattering in the s± superconductor is negli- 
gible. In this case bulk superconductivity is described 
by quasiclassical Usade l eq uations^ with boundary con- 
ditions derived in Ref. 123 . The conventional proximity 
effects were extensively explored within this approach in 
Ref. This description was later generalized to conven- 
tional two-band superconductors in Ref. [3ll This model 
have been already used to describe some anomalous prop- 
erties of the s/s± interface in Refs. l22l 

Both superconductors are described by the gap param- 
eters, A(x), and the impurity averaged Green's functions, 
which have regular and anomalous components, G(x, w) 
and F(x,lj), with G 2 + |F| 2 = 1, where u = 27rT(n + l/2) 
are the Matsubara frequencies. In the following, we 
will use subscript "s" for the s-wave superconductor 
and subscript "a" for the a-band of the s± supercon- 
ductor and skip subscripts in relations applicable for 
both superconductors. Further, we employ so-called 
parametrizatio n 30 ' 31 in which the function $ = uiF/G 
is used instead of F. In this case G = ujj y/cu 2 + |$| 2 . 
For the s-wave superconductor (— d s < x < 0), the equa- 
tions for the Green's functions G s and $ s and the self- 



consistency equation are: 
D s d 



2uiG s dx 



2.TJ2 



-A s , (4a) 



A, 



A s ln^=0 (4b) 



Correspondingly, for the s±-superconductor, < x < d± 
we have 



2ujG a dx 



G 2 d<f> a 



a dx 



2ttT > X a g — — 



= A c 



(5a) 
(5b) 



where \ a p is the coupling-constants matrix and a, j3 are 
the band indices. In Eq. (|5"aj) we neglected the interband 
impurity scattering. For the case of s± superconduc- 
tor we consider here A1A2 < 0. This is realized when 
A12, A21 < 0. The diffusion coefficients -D{ s ,a} are related 
to the conductivities ct{ s ,q} as <J{ s , a } = e V{s,a}Ds Sj0l \, 
where V{ s , a } are the normal densities of states (DoS). 
The ratio of the off-diagonal coupling constants is given 
by the ratio of partial normal DoSs, A Q| g/A^ a = v^fv a . 
It is convenient to normalize all energy parameters (ui 
and gaps on both sides) to the same scale irT c . We 
also introduce the coherence lengths £ Q = \J D a /2irT c 
and £* = 'D s /2ttT c (note that £* is related to the 
bulk coherence length of the s-wave superconductor by 

The bulk equations have to be supplemented with the 
boundary conditions at the interface separating two su- 
perconductors. These conditions relate the Green's func- 
tions and their derivatives at the interface and can be 
written as 2 ^* 



£a Go. 



d<5> s 
dx 

d®g 

dx 



X)^(*a-*.), (6a) 
--— ^($ S -$ Q ), (6b) 

iBa 



for x — 0, where a is the band index. Here the coupling 
parameters, 73 a and 73 a are proportional to the partial 
boundary resistances Rbu, 



IBa = 



RBa 

ps£' 



RBa 
Pa£,a 



(7) 



where P{ s , a } — l/f{s,a} are the bulk resistivities. We will 
also use the ratios of these parameters 



7a 



iBc 
iBc 



Pa£,o 
Psts 



(8) 



that are bulk parameters characterizing the relative 
"metalicity" of the s-wave superconductor and a band. 
In particular, large 7^ implies that the s-wave material 
is more metallic than the a band on the s± side. The 



4 



parameters j a , \ a p, and £ Q are not fully independent, 
the ratio of 7 a obeys the following relation 



7i 

72 



^2^2 -^12^2 



^l£l <^2l£l 



(9) 



The conditions at the external boundaries are &' s (— d s ) = 
and $' a (ef±) = 0. 

The supercurrent flowing through the interface be- 
tween the s-wave superconductor and a band is deter- 
mined by the Green's functions ${ SiQ } at the interface 
as 



Im[*:$ a ] 



u>0 



Vw 2 + |$ s |V w2 + l*al 



(10) 



with ^4o = V( e Ps£s)- Substitution of the zero-order ap- 
proximation for the Green's functions, ${ s ,q} ~ A/ S)C( jq, 
gives the well-known Ambegaokar-Baratoff result^ 2 - for 
the partial Josephson currents proportional to sin <j> with 
different signs corresponding to the signs of A q0 . Here 
and below we assume for definiteness that <j> is the phase 
shift between A s o and Aio. To find the sin(2^>) term 
in the Josephson current one has to go beyond the 
zero-order approximation and evaluate corrections to the 
Green's functions due to the interface. We discuss these 
corrections in the next section. 



ary resistance Rsa by the standard expression 

A s o|Aq,o| 



Ej,a — 



2e 2 R% ^ 



% V^+a^^ 2 +|a q0 | 



(ii) 



This also means that the total Josephson current flowing 
through the boundary is always zero in the ground state. 
The corrections can be presented in the form of Fourier 



lm[A] 




FIG. 2. Illustration of the bulk gap parameters, A 3 o and Acq, 
and interface-induced corrections to the average s± gaps, A Q , 
for a general TRSB state. The gap parameters are presented 
as vectors in the complex plane. Decompositions of A s o and 
A a into the amplitude and phase components are also illus- 
trated. 



IV. PROXIMITY CORRECTIONS IN THE 
WEAK-COUPLING LIMIT 

In the case of weak coupling between the s and s± 
superconductors, jBa 3> 1, the contact-induced correc- 
tions to the gaps and Green's function can be treated 
as small perturbations, A^ s ^ a y(x) = A{ sa } + A{ sa y(x), 

®{s, a }(x) = &{ s ,a}o + ®{s,a}{x)- As a zero-order ap- 
proximation, we consider a general complex case with a 
finite phase difference <p between the bulk gap param- 
eters Aio and A s0 , see Fig. [2j For the aligned states 
such perturbative calculation has been reported in Ref. 
I22I . Without loss of generality, we assume A s q to be real. 
The small corrections $/ S)Q ,i(o;, x) and Ai St0l \(x) can be 
computed analytically in the linear order with respect 
to 1/jBa- Similar calculation for several types of junc- 
tions using somewhat different approach has been done 
in Ref. |H. The details of these derivations are described 
in Appendix [3] In the TRSB state the solution for cor- 
rections exists only if the partial Josephson energies in 
the linear approximation exactly compensate each other, 
Ejx = Ej2- Here Ej a are related to the gaps and bound- 



expansions. For the s-wave superconductor <f> s (ui,x) — 

X)m=0 $ s,m(w) COSfc m IE, A s (x) = J2m=Q ^s,m COS k m X 

with k m — m/K/d s . The Fourier components of the 
Green's functions computed in Appendix IA II are given 
by 

- A s , m (2-S m )e s ,J{dsC a ) 

x y^ V^+^o A a0 -A s0 

where ^ = £ 2 , A A s0 /^ 2 + A 2 , £ 2 A = D s /{2A s0 ), and 
S m = 1(0) for m = 0(m > 0). Here the first (bosonic) 
term is induced by the correction to the gap parameter 
and the second (fermionic) term is the direct response 
to the boundary perturbation. In the complex state the 
responses of the gap parameter are different in the ampli- 
tude and phase channels. As A s q is selected real, these 
channels correspond to the real and imaginary parts of 
the gap correction, A s = Af + iAj, . The Fourier compo- 
nents of A^ and Af computed in Appendix IA II can be 
presented as 
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^ = 2-kT w» (2-3 m )ff A /(rf a &) Re[A a0 ]-A s0 



1 | A 2 uj 2 {nm£ StA /d s 



WithZ -" 27rT S(- 2 + A ?0 )3 /2 ( v A -- v /- 2 + A^ /A s0 + (^aM)V 
A/ = 2^T 2^ A /(4g.) 1 Im[A a0 ] 

with Z^ m = 2vrT 



f v 7 ^ 2 + A 2 ^/w 2 + A 2 /A s0 + (nm£ sA /d s y 



For the s± superconductor the corresponding 
Fourier series are $ Q = Ylm=o & a ,m cos q m x, A Q = 
ELo A Q , m cosg m a; with g m = mn/d±. Detailed 
derivations of the Fourier components are presented in 
Appendix I A 21 and the result for $ Qjrra can be written as, 



1 ■ 



£2 „2 



(14) 



where 



^=£ 2 ,A|A a o|/v^ 2 + |A Q o| 2 , £ A = D a /(2\A a0 \). 

Here, as in Eq. (I12[) . the first term is induced by the 
correction to the gap parameter and the second term is 
directly induced by the interface. For a general com- 
plex state, the corrections have to be split into the 
amplitude (along A Q o) and phase channels, Q a ,b.m — 



Oi,b. m 



(2 - j m ) eiL/ (d±£ a ) 

IB a 1 + £L,?m 



v/c^ + IA^ol 2 /A° n -A tt0 



A 2 



(15) 



with the matrices U m ^ a p 



[Wa/3 



ya 

a , m 



where 



' L( w 2 + |A Q o| 2 ) 3/2 (l+£ 2 ,^) w . 



T 



T 



V ^+|A^f (1 + e 2 2J w 



= A Q( ^ — A 1 5 a p and A is the largest eigenvalue of the 
matrix \ a f). In the two-band case the explicit formulas 
for w a f3 and U^ a/3 are given in the Appendix IA 21 Eqs. 
(|A"T9)) and l|A"27|i . 

In summary, Eqs. (fT2"]) and (fT5|) give corrections to the 
gap parameters and Green's function for the s-wave su- 
perconductor while Eqs. (fl4|) . (fl"5j) . and (|T6|) give corre- 
sponding results for the s± superconductor. These cor- 
rections will allow us to derive in the next section a gen- 
eral result for the second harmonic of the Josephson cur- 
rent that determines the width of the TRSB region. 



SECOND HARMONIC OF THE JOSEPHSON 
CURRENT 



where A" (A* ) is the projection of A s0 along A q0 (into 

perpendicular direction), as illustrated in Fig. [5J The 
gap corrections A£>£, are related to $^ m as 

K, m = 27TT Y: U^ ap "^A '^3/2 ' ( 16a ) 
0,a»O (^ 2 + |A/3o| 2 ) 

A* m = 2^T £ <^ > f^'T , a , (16b) 
/3^o V^ 2 + I A/3o 1 2 



We already mentioned that the linear order with re- 
spect to the coupling strength oc is not sufficient 
to determine the range of parameters where the TRSB 
state is realized. As discussed in Sec. lU in the weak- 
coupling regime this range is determined by the term 
cx sin(2</>) in the Josephson current that appears only in 
the quadratic order. In this section we derive microscopic 
expression for this term using corrections to the Green's 
functions presented in the previous section. 

For arbitrary coupling the current density flowing 
through the interface between the s-wave superconduc- 
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tor and a-band is given by Eq. (|10[) and in ground state 
3 = E- 7 " = Aal = 0. 

a 

Here the parameter 1 has dimensionality of energy. To 
find the second-order term in j, we have to expand the 
right hand side of Eq. (|10[) with respect to small correc- 
tions to <E> S and <I> Q . This gives I « 1^ + Z (2) where the 



term 



2T« = 2nT 



Im[A* A Q0 ] 



a j^0 jBaV^ 2 +\^ao\W^ + \ A so\ 2 



(18) 



corresponds to the standard main-order Josephson cur- 
rent and 



1^ = 2ttT 

a,u>>0 



Im[$;A a0 + A; $ a ]-Im[A; A a0 ] + 
7saVw 2 + |A Q o|V w2 + l A «o| 2 



(19) 



is the second-order term which is determined by the lin- 
ear corrections to the Green's functions due to the in- 
terface perturbations, ${ S)Q }, considered in the previous 
section. Using these results, we can present the correc- 
tions at x = in the form 



#*(0)=V;j-« A, + (-l)"|A a0 |cos^ 



iBa 

(-l) Q |A Q0 |sin< 



(20a) 



(20b) 



a/3 



IBP 







a/3 „, 
7B/3 



(20c) 



where the response functions of the s-wave superconduc- 
tor in the amplitude and phase channels, ^"^(w), can 
be explicitly written as 



Ta = ewe ^ 2 + A 2 f; 2 -a, 
• - - - ' i 



tanh(d s /£ S;W ) ,/u; 2 + |A q0 | 2 



2ttT x 

x £ 



1=0 1 ^cd'-m 



% (uf + A 2 ) vM + IAool 2 1 + £W*4 ' 

(21a) 

WC y/^ 2 + A 2 ^ 2 



tanh (4/6,a;) ^/w 2 + |A q0 | 



m=l " ' * s > u; m 

x ML v i ^WC^Cj f 2 ib) 

x ~ ^? + ia q0 i 2 i + e a<m ki' 



with defined in Eq. (fT5)> . The response functions of 



the s± superconductor, ^'^(oj), are given by 



« ^ 2 + |A a0 |% ~ 2-6 

J-„r = - , ; , r , ; = d a p+ y 

sO m=0 



tanh(d±/6W v^ 2 + A 2 



£2 „2 



wi>0 



/3 ( W 2 + |A /3 o| 2 )v / ^A^ 1+^^^ ' 

(22a) 



Uu/tc, V" 2 + |A 



tanh (d±/£ a ,u) yjuj 2 + A 2 



2 ^ 1 

1=1 T ^* Vm Wl >0 



(22b) 



where the matrices a/3 are defined in Eq. (|A27|i . and 
the term 



-p'i- 



2ttT 



£ 



2d ± w 12 \A 20 \ ^ ^ 2 + A 2 



(23) 



describes the contribution from the uniform phase correc- 
tions A^ given by Eq. (|A29[) . The different quantities 
entering Eq. (|19[) can now be expressed as 

Im[A^A Q o] = -A a o|A a o|(-l) a sin^ 

Im[$:A Q0 ] = |A Q0 | [-^%A sO (-l) a sin0 

p Ibp 

- (jVj-T*^ |A^ o |(-l) a+/3 sin0cos0 
Im[A: $J = A s0 £— [(-l^lA^olsin^ 

A 7B/3 



^a/3 ~ J A s0 SUl (j) COS ( 

Re[A: ^ = A s ^^ A - + ^^l A ^' C ° S ' 



7S/3 



^=-(-1)1^^ co S ^(-l)^|A, | 



7S/3 
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We can see that I 1 - 2 ' contains terms proportional to 
sin</> and sin^cos^ = |sin(20), l^— sin + 
J^ 2 - 1 cos ^ sine/). The terms oc sin</> just slightly shift lo- 
cation of the transition line. The transition order and 
possible width of the TRSB region are determined by 
the terms oc sin <fi cos <fi. Collecting such terms in I, we 
obtain 



\A a0 \A s0 K a (cu) 



n a {io) = J2 



(24) 



s0 U 2 + A 2 ^ 



L 
A s0 



IB A 



/ . 2 T^a 



7 B/3 |A q0 | \ + |A a0 | 



J" 



a/3 



The term corresponding to the uniform phase tilt T^ , 
Eq. (j23|) has special meaning and it is useful to evaluate 
it explicitly, 



7 (2) _ 

u 4>,0 ~~ 



2nT 



E 



A s0 |Ai 



7i^,a|Aiq|/6 
d±wi2|Aio||A 2 o| 



2e 2 
= — 



Ps£,s 



7,1 



d±viWi2\A w \\A 2 Q\ 



(25) 



where we used the relations 7i£ 2 A | Aio|/£i = 
hp\D 1 /{2p s ^ s ) = h/{2e 2 V\p s £ >s ). This gives the follow- 
ing contribution to the second harmonic of the Josephson 
current 



.(2) 

3<l>,o 



-JJ,i 



2d ± i/iwi 2 |Ai ||A 



2d 



(26) 



The quantity £\2 = d±^iWi2|Aio||A2o| in the denomina- 
tor represents the interband coupling energy. Therefore 
this result exactly corresponds to the result of the frus- 
trated Josephson junction model. This term, however, 
dominates only in the case of small interband coupling 
energy, when the parameter w\% is very small. 

To gain a further insight on the structure of the second- 
harmonic amplitude J^^ 2 \ we present it explicitly as a 
sum of terms corresponding to contributions from the 

corrections to the s-wave and s± Green's functions COm- 
^l 

ing directly from the interface ( J7" b J ) and via gap param- 
eters (J™). 

M2) _ T (2) ^(2) _^ ^(2) 
U Uh a ~t~ u a s 



J, 



fc,pm 



J 



(2) 



A,pm 



T (2) 



(27) 



(2) 

We already considered above the last term, , that is 

(2) 

part of J7 A p m coming the uniform phase response of s± 
superconductor. This term requires separate treatment 
leading to Eq. ([25]) which corresponds to the result of the 
frustrated Josephson junction model. As for the other 
terms, the s-wavc components arc given by the following 
explicit formulas, 



7(2) 

Jb,s 



7(2) 



- 27rT E E~ 



A Q o 



A 2 



u>0 



>.S' 



lB«^u 2 + |A a0 |2 J W 2 + A 2 tanh(d s /^)' 

y2 00 / y2 y2 



J s,0 



Y^ m = 2nT J2 



7 4> go, 

A a Q 



(28a) 
(28b) 



,a,u»0 iBa y/w 2 + lAaOp [y/ U) 2 + A'j /A s0 + Ca.A) 



Y a , m = 2irT J2 



A a0 u> 2 



rn , a , a; > ^jBa V^ + |A Q0 | 2 ^ 2 + A 2 ) (^ 2 + A 2 /A 
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while the s± components can be written as 

a 2 |a q0 | 2 



J, 



(2) 



6.pm 



7 (2) 
'-'A.pm 



-2ttt y ^ 



7B a 7sa (w 2 + A 2 s0 ) (u 2 + |A q0 | 2 ) tanh(d±/£ a , a 



7rT r 



-A 



oo 

- X O, a U O,<x0l0t0 X O,0 + 2 ^2 ( X ™,a U t, a /37pt,P X t,p ~ X m,aUm,<x0l0£.0 X m,c?) 
rn — 1 

l 



X' 1 



2ttT x ^ 



(29a) 
(29b) 



o V^ 2 + A 2 V^ 2 + lAaol 2 + |A Q0 |a a ql 



2ttT 



fo v/^ 2 + A 2 ( W a + |A Q0 | 2 ) ( V ^+]A^+ |A Q0 |el,X) 



We note also a useful relation for the combination 7 l a£ l a 
entering Eq. (|29bl) . 73^3 = t;* s v s /vp. In summary, Eqs. 
(|23|) . (|2"5|). and (f2T))) give general expressions for the com- 
ponents contributing to the second harmonic of Joseph- 
son current in Eq. P7|) . Even though these formulas are 
rather cumbersome, they are suitable for numerical eval- 
uation of the second-harmonic amplitude for arbitrary 
parameters of superconductors and interface. In the next 
section we analyze these terms for practically important 
particular case in which much simpler analytical results 
can be derived. 



A. Analysis of terms for low temperatures in the 

case d 3 < £ s and A s <C |A a | 

Unfortunately, general formulas derived in the previous 
section are rather cumbersome. To understand better the 
relation between different terms and their absolute val- 
ues, in this section we evaluate them at low temperature 
and for the most interesting case of weaker s-wave su- 
perconductor, small d s and large d±. In these limits it is 
possible to derive simple analytical results for the most 
important terms. 



1. Terms J$ and J£ 

For very thin s-wave superconductor the dominating 

(2) 

in 1 / d s order term is coming from J and the m = 

(2) 

term in the amplitude part of ' which we will notate 
as T (2) 




A Q o 



A 2 

^sO 



- 7BqV /^ 2 +|a q oIV (^ 2 +a 2 o) 3 / 2: 



a q0 a 2 



.10 



7W^ 2 +|Aao| 2 (w 2 +A 2 ) 



3/2 



r 



where in the last formula we used a compensation condi- 
tion for the Josephson energy near the transition point. 
In the limit A s0 <C |A a0 | integrals converge at ui ~ 
A s o and one may think than it is possible to replace 

|A a0 | 
% /w 2 +|A q0 | 2 

as / °° {u]2+/ ^ o)3/2 = 
two terms in the sum 



1 under the frequency integrals. However, 

2 

1, within this approximation the 



7 (2) 



j(2) 



J 



(2) 



A,s,0 



exactly compensate each other, 




and therefore they must be evaluated in higher order with 
respect to A s0 /|A Q o|- To proceed, we introduce the def- 
initions 



(-l) a |A a0 | 



- jBay/u 2 + |A a0 | 2 

with L = E Q ( _1 ) Q ' llBa and 



L + Li(w), 



(-1)° 



lAool 



V^ + IA^oT 



which allows us to represent 



J £=-j-J dujL ' Z{uj 

1 / poo 

and rewrite Jq 2 ] as 



1 



A 2 



(^ 2 +A 2 ) 3 / 



( W 2 



A|o 

+" A 2 ) 3 ' /2 / 



ft 



(2) 




duLi(w)- 



A?) 3/2 



duoL\(uj) 



A 



(lu 2 + A 2 ) 3 / 2 
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The dominating contribution is coming from the second 
term and evaluating integral, we finally obtain 



7- (2) ~ —a 



A 



2 In 2 - 1 



(-ir 



In 



(i+fcl) 



7Ba|A Q | J 



2 In 2 
|Ai||A 2 | 



7bi7S2 



(30) 



It is interesting to note that this l/d s term only exists in 
the asymmetric case, it vanishes for identical s± bands. 



2. Term 3$^ 



The term J^ m in Eq. ([29a] for T = and d± > £ 
becomes 



7 (2) . 



A^ |A Q0 | 



IBolIBol 



2^/4- 



(31) 



For A s o <C |A Q o| we can evaluate the frequency integral 
leading to the quite simple result 



7 ( 2 ) 



7T X - A^q 

2 ^ jBajBa 



±a0 



(32) 



In most cases this is actually a dominating term which 
may be used for an approximate evaluation of the to- 
tal second-harmonic amplitude. It has only linear order 
with respect to A s o, while other terms are proportional 
to A 2 . Also, for a typical contact j a = jBa/lBa ^ 1 
due to semimetalic nature of iron-based superconductors, 
which enhances J*"~p m terms in comparison with jj 2 ] 

(2) 

terms. The negative sign of J b ; pm implies the continuous- 
transition scenario and the existence of the TRSB state. 

(2) f2) 

In particular, comparing J b J m with term JqJ, Eq. (|3"U)) . 
we obtain, up to a dimensionless function of the ratios 
iBillBi and |Ai |/|A 20 |, 



7"W 



d s |A 



7i 



101 



£ A s0 ' 



<r(2) 



exceeds 



which means that the s-wave term 

(2) 

J b , m only for an extremely thin s-wave layer d s < 
£*A s0 / (7i|Aio|). Another factor further enhancing this 

(21 

ratio is that different s± bands contribute to J h ) ' with 



the same sign while their contributions to __7 
compensate one another. 



(2) 



6,pm 

partially 



3. Termji 2 ^ 

Finally, we obtain the limiting form for the most com- 
plicated term in Eq. ([29b| . The quantities S^f Q 

that determine the matrices U^ a g in Eq. (|A27j) can be 
evaluated at T = exactly: 



^(1/2), 




(33a) 
(33b) 
(33c) 



where a m — {^rn^ a /^/d±) and ip(x) is the digamma 
function. The quantities X^ a can be evaluated approx- 
imately in the limit A s q *C |A q q|, 



1 



1 



In- 



4IA 



n0 



7Sq|A q0 | l + a m V A s0 
1 ln(l + a m ) 



l-a„ 



In 



1+a, 



(34a) 
(34b) 



7Ba|A a o| a m 
Assuming that the values for £ Q , |A Q o|, and 7 Q for diffcr- 

can be roughly estimated 



ent a are close, the term J L 
as 



(2) 

A,pm 



A.pm 



TrT r 



A 2 

^s0 



7si7BiA 2 



(35) 
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and we can see that this term is typically smaller 

(2) _____ 

than J~ bpm , Eq. (f3"2"l) . by the ratio A s0 /Ai (assuming 
Ai /7rTe~l). 



B. Region near T c s for T<f < T c 

Near the transition temperature of the s-wave super- 
conductor all terms contributing to decrease as A 2 . 

(2) 

In particular, the dominating term Sf b ' behaves as 



7 (2) 




A a o| jBalBa 



forT^T c s -0. (36) 



This behavior has an important consequence: the width 
of the TRSB state shrinks with increasing temperature. 
On the other hand, the weak coupling approach breaks 
down when the temperature is too close to T£ when the 
correction to the s-wave gap becomes comparable with 
its bulk value. 



VI. WIDTH OF TRSB REGION 

To analyze the width of the TRSB region in the weak- 
coupling regime, we represent the supercurrent flow- 
ing through the interface in the form j = ^ a j a = 



10 



Ao (J^ sin 4> + cos 4> sin . The transition roughly 
corresponds to the vanishing of the first harmonic 
which we can represent as 



iBa 



(37) 



with 



,- r i I AgO | Asp 



= 2ttT 



w>0 



V^ + IA^IV^ + IA^ol 2 ' 



We remind that the TRSB state only exists if J {2) < 0. 
In this case, which is realized for our system, we can write 
the condition for the TRSB state range as 



fj,i 

TBI 



fjg 

7B2 



< 



\JV>\ 
A s0 



(38) 



This formula together with microscopic results for 
of the previous section represent the main results of this 
paper. For fixed jg} the transition from aligned to TRSB 
state occurs at the following values of 7^2 



1 

7-B2 

As scales as -■' J 



± 



fjalBi A s0 /j,2 

7 B ", The width of the TRSB region 
can be conveniently characterized by the parameter 



l 2 B2^lB2 ~ A 7 S 2 = 



27| 2 J (2) 

A s o/j,2 



(39) 



which depends only on bulk properties of the superconduc- 
tors and does not depend on the boundary resistances. 34 
In particular, at low temperatures and for A s q <C IA20I, 
using the asymptotic fj j0l « ln(4| A a0 |/A s o) and keeping 
only the dominating term in J^ 2 \ Eq. (|32|) . we obtain 
the following estimate 



A7. 



B2 



7T In 



4|A 2 



■E 



7« 



(40) 



We can see that this width only weakly depends on the 
value of the s-wave gap and is mostly determined by the 
parameters j a . Using the definition of jBa and 7a , Eqs. 
([7|1 and (|5]), we immediately obtain a simple estimate 
for the spread of the partial boundary resistance ARb2 
within which the TRSB state exists, 



AR 



B2 



it ln- 



4|A, 



A s0 



Pa^a 



I A q o| ' 



V(ln^i 
which is mainly determined by the products p a £, a 



(41) 





X 12 =-tJ.4, X2,= 


■° y ' / 


0.02 


"71= 4, d s =£* 


/ ft* 


0.04 




/■ /■' 


0.06 




ft 


o.ox 












0.10 




// (a) ; 


0.12 






0.14 




. • • r <>.?t. 






• TJ=0.5T C 


0.16 







0.2 0.3 

T/T„ 







: (b) \ 




T^=0.3T C 




T^=0.5T C 


I.. A 


b.o 0.1 0.2 


0.3 0.4 0.5 



T7T„ 



FIG. 3. (a)Temperature dependences of the second-harmonic 
amplitude for two values of the ratio /T c , 0.3 and 0.5. 
Dashed lines show the term j£ m only. Other parameters 
are shown in the plot. We also assume An = A22 = 0. (b) 
Corresponding temperature dependences of the TRSB width, 
Eq. 



Figure[3]illustrates the temperature dependences of the 
second-harmonic amplitude and the width of the 

TRSB region for the representative parameters listed in 
the left figure and for the two values of the ratio T c s /T c , 
0.3 and 0.5. In Fig. Eta) we show for comparison both 

the full amplitude of the second harmonic and the term 

(2) 

i7j pm only (dashed lines). We can see that this term 

typically accounts for 80% — 85% of the total amplitude. 

f 2) 

The remaining part mostly comes from the terms J7/ 



(2) 

and J\ J . For used representative parameters the contri- 

(2) 

butions from the s-wave terms J7*\ s are negligible. We 
emphasize the shrinking of the TRSB width illustrated 
in Fig. 13(b), A 7B2 oc y/T* - T near T c s . This means 
that in some range of parameters, the transition from 
the aligned to TRSB state may be observed as a func- 
tion of temperature, as in the case of the 0-n transition 
in SFS junctions^ 



A,pm 



VII. DENSITY OF STATES OF S-WAVE 
SUPERCONDUCTOR WITHIN THE TRSB 
REGION 



Contact with s± superconductor induces specific fea- 
tures in the s-wave density of states (DOS). For the 
aligned states in the case of thin s-wave layer, the s± gaps 
aligned with A s0 generate positive corrections to DoS, 
while anti-aligned gaps generate negative corrections? 2 ^ 
The latter negative features can be used to identify s± 
state. 

In this short section we consider evolution of contact- 
induced features of the s-wave DoS across the TRSB re- 
gion. To find the density of states, we have to perform 
analytical continuation of the Green's functions to real 
energies iu) — s- E + i8. The normalized DoS is related to 
the real-energy Green's function by the standard expres- 



11 




FIG. 4. Evolution of the proximity-induced features in the 
s-wave DoS within the TRSB region with increasing angle 
cj> between A s o and Aio. Curves are vertically displaced for 
clarity. Dashed lines show the bulk DoS. 



N(E, x) = Re 



E 



sjE 2 - $(E,x)&(E, 



(42) 



where &(E,x) = $*(—E,x). Expanding § s (E,x) and 
taking into account that A s o is selected real, we obtain 
the proximity-induced correction to the s-wave DoS 



5N s (E,x) fa Re 



EA s0 (*l(E,x)+* B (E,x) 



2(E 2 



A 2 ) 3/2 

^sOJ 



(43) 



obtain for these contributions at real energies, 



t 



(2-6 m )e s , A /(dsO 



v / E 2 -A 2 s0 /A s0 + i (nmi sA /d s 



E 



V E2 - A lo Re[A Q0 ]-A s0 



- ^\A a0 \ 2 -E 2 



JE 2 -A 2 A R 

V sO s.n 



v/^ 2 -A 2 + iA s0 (nm^ A /d s ) 2 ' 



where A R m is given by Eq. (fUSa ). These results together 
with Eq. (|43|) determine the shape of the DoS correction 
for arbitrary parameters of superconductors in the lin- 
ear approximation with respect to the coupling strength 

1/jBa- 

For the important case of small d s we can keep only 
the uniform m = term in the Fourier expansions leading 
to a simple result similar to one for the aligned stated 2 ., 



^ Re[A a0 ]- 



Vs() 



01 



£ 2 



(44) 



with 



^ 2 

d f* 



2e 2 v s d s R% ' 



Therefore, the correction to the s-wave DoS in the case 
of thin s-wave layer is given by, 



SN s (E,x) 



EA s0 



(E 2 



v2 



x Re[Aao] - As0 Q(\^o\-E), (45) 



VIAaol 



-E 2 



where Q(x) is the step function. If, as before, we de- 
fine that phase shift between A s o and Aio as </> then 
Re[Ai ] = |Ai O |cos0 and Re[A 2 o] = — |A 2 o|cos0. 

Figure @] illustrates the evolution of this correction with 
increasing angle <f> for representative parameters. Two 
limiting cases cj> = and 7r correspond to aligned states 
in which the aligned and anti-aligned gaps induce asym- 
metric peak and dip correspondingly. With increasing 
phase, the peak smoothly transforms into a dip and vice 
versa. In the maximally frustrated state for <p — n/2 the 
DoS has two small dips. Note that the DoS correction is 
obtained with the linear approximation with respect to 
the coupling between the superconductors, which some- 
what overestimates the amplitude and sharpness of the 
peaks j^ 2 . 



The correction to the Green's function can be 
represented again as a Fourier series, <& S (E, x) — 
Sm=o ®s,m(E) cos (rmrx/d s ). In the Matsubara presen- 
tation, two contributions to the Fourier components <fr s . m 
are given by Eqs. (|A8I) and (|A9I) . Using these results, we 



VIII. SUMMARY AND DISCUSSION 

In summary, we evaluated the range of parameters 
where the TRSB state is realized for the interface be- 
tween s-wave and s± superconductors using the simple 
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microscopic theory. This state appears when the partial 
Josephson energies almost completely compensate each 
other. The width of the TRSB region is determined by 
the sin(20) term in the Josephson current, which appears 
in the second order with respect to the coupling strength 
between the superconductors. We found that the dom- 
inating contribution to this term is determined by the 
direct boundary correction to the Green's function of the 
s± superconductor. This term is missed by phenomeno- 
logical models of the junction. The width of the TRSB 
region shrinks with increasing temperature giving possi- 
bility to detect the transition from the aligned to TRSB 
state as function of temperature. 

The main purpose of this paper is to establish factors 
which determine the width of the TRSB region at the 
s/s± interface in the simplest possible situation acces- 
sible for full analytical analysis. Even in this relatively 
simple case the analysis occurred to be very nontrivial. 

Several factors may have quantitative influence on re- 
sults reported in this paper and complicate their appli- 
cations to real iron-based superconductors: 

• Most of these materials have more than two bands 
(up to five). This is not a crucial complication. Our 
consideration can be directly generalized to arbi- 
trary number of bands. 
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Appendix A: Derivation of boundary-induced 
corrections in the weak-coupling limit 

In this Appendix we derive corrections to the Green's 
functions and gaps in the linear order with respect to 
the coupling parameters 1/jBa- As a zero approxima- 
tion, we consider a general complex case when there is 
a finite phase difference 4> between the bulk zero-order 
gap parameters Aio and A s o as shown in Fig. [2] Corre- 
spondingly, the phase difference between A20 and A s o is 
cj) — TT. This calculation covers both aligned states when <f> 
equals or 7r and complex TRSB states with < cf> < ir. 



• Due to the very short coherence length, most 
iron-based superconductors are probably in the 
clean limit. In the cleanest materials a significant 
anisotropy of the gap2& and even the gap nodes^l 
were revealed. On the other hand, in several other 
compounds isotropic gaps within the bands were 
found 38 , which probably indicates substantial in- 
traband scattering. The presence of a significant 
gap anisotropy and nodes does not contradict an 
overall picture of the s± state because what mat- 
ters most is the average gap inside the band. This 
means that the TRSB state is also expected within 
some range of parameters at the interface between 
conventional and clean s± superconductor. How- 
ever, our calculation of the width of this region is 
not directly applicable to this case. 

• In compounds with strong impurity scattering one 
can expect some interband scattering, which was 
neglected in our model. As this scattering sup- 
presses s± state, it can not be too strong. The 
main effects of this scattering are suppression of 
the s± gap parameters and emergence of the sub- 
gap states. These effect may have some influence 
on the location and width of the TRSB region. 

An accurate description of these factors requires special 
considerations that will further complicate the theoret- 
ical model. Nevertheless, we expect that a qualitative 
picture of the transitional region will hold within a more 
realistic framework. 



1. s-wave gap and Green's function 

We start with calculation of corrections to the s-wave 
Green's functions and gap parameter, <I> S and A s . From 
Eq. (|4a|) we obtain that the first-order corrections to the 
s-wave Green's functions obeys the following equations 

C^-$ s = -A s , (Al) 

where & >u = D./ (2^+^) = £ s 2 >A A s0 /V^+A^ 

and £ 2 A — D s /(2A s0 ). Without loss of generality, the 
zero-order gap parameters A s o can be selected real. In 
this case the self-consistency condition for the linear cor- 
rections can be written as 

In the boundary condition for d$ s /dx, Eq. (j5a]), we can 
neglect in the right hand side differences between $'s and 
A's and approximate A's by their bulk values. This gives 

V^+Aj da ^ jBa ^ 2 + |A Q0 | 2 

at x = 0. Note that, in general, A s , $ s , and A Q o are 
complex, A q0 = | A q0 | exp(i0 a ), fa = fa fa = <t> ~ T- 
They are real only for the aligned state. 
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To solve Eqs. (|A1I) and (|A3|) it is convenient to split 
$ s into the two contributions, <t s = <b s b + $ s .a, where 
$ Si 6 is induced by the boundary condition and $ S! a is 
induced by the gap correction. The first contribution 
$ s b can be found from the following equation and the 
boundary condition 



£A - = o, 



(A4a) 

g^, b (0)=E^- rSP% ( A *o-A s0 ), (A4b) 
while the second contribution, $ S ,A, obeys 



(A5a) 
(A5b) 



The solution § S} i,(x) of the linear equation (|A4aP with 
the boundary condition $^ b = at a; = — d s is given by 



^ s .b(x) = C s . 6 cosh 



i + d s 



(A6) 



where the constant C S) & can be found from the boundary 
condition at x = 0, Eq. (|A4b|) . 



a 



s,b 



SWh(d s /£ s , u ) ^ jBa^UJ 2 + \A a0 \ 2 

leading to the following result 

£ SiU cosh [(x + d s ) /£ s , w ] 



(A Q o — A s q) 



$a,f>(z) 



£*sinh(d s /£ SitJ ) 
x y_V^^2=(A a0 -A a o). (A7) 



We compute $ s ,a and A s using the Fourier ex- 
pansions, $ s ,a( x ) = Em=o ®s,A,m cos fc m a:, A s (a?) = 
J2m=o As,m cos fc m £ with fc m = mir/d s . The Fourier 
components of $5^(2;) can be computed explic- 

itly from Eq. ([A?]) . 



_ (2-^m)^/K^) y^ 2 +A 2 A a0 -A, 

1 + 4" V w2 + |Aao| 2 7Ba 

(A8) 

Eq. (|A5aP immediately gives the following relation be- 
tween the Fourier components $ s ,A,m and A SjTO 



$s,A,m = 



A. 



1 i C2 1.2 ' 

- 1 - ' Ss.u^ra 



(A9) 



Substituting this result into the self-consistency condi- 
tion (IA2[) and splitting it into the real and imaginary 
parts, we relate A SjTO = Af + iAf to l> s , 6, m = 



2^E 



^2 i , ,2 «s,oj fc 



1 s0" 



-CJ 



2rh-R 



w> ' (^ + A2 ) 3 / 2 

£ 2 k 



A£ m = 27rT£ 



^ ( W 2+A2 ) 3 / 2 ' 

(AlOa) 



i f 2 k 2 4> / , 

2tfV^= ^ ?, Af m = 2^rV s ' h ' m 

^ v^ 2 +a 2 i+e.,fc 2 s ' m ^ 



s.uj m 



.>oV^+A^ 
(AlOb) 



The components Af m and Af m describe responses of 
the s-wave gap parameter to the interface perturbation in 
the amplitude and phase channels. Eqs. (|A8|) . (|A9j) . and 
()A10|) already provide a formal solution of the problem. 
As the left hand side of Eq. (|Al0b[) vanishes for m = 0, 
solution for the imaginary part exists only if 

y .Ski =o 

giving the condition 

\ - 1 Im[A a o] 



0. (All) 



Uncertainty in Af reflects the phase-rotation invariancc 

and we can select Af = 0. As the Josephson energy 
between the s-wave superconductor and a band, Ej a is 
given by Eq. ([TT]) and A Q o = |A q q| exp(i(/) a ), we can see 
that the condition (| Al 1|) simply means 



0. 



Since the partial Josephson currents are proportional to 
Ej a , this condition implies that the total Josephson cur- 
rent flowing through the boundary is always zero in the 
ground state. For two bands the condition for realization 
of the TRSB state in the linear order with respect to the 
interface transparency is simply Ej\ = Ej% . To establish 
an accurate range of parameters within which the TRSB 
state is stable, one has to go beyond the linear order. 

Using the expansion (|A8[) . we obtain the explicit pre- 
sentations for Af m and Af m given by Eqs. (fl"3"|) of the 
main text which in turn determine the term $ s ,A,m of the 
Green's function, Eq. (IA9[) . Since we already derived the 
result for $ s ,b,mi Eq. <1A7[) . we now have all corrections. 

The full analytical formulas are somewhat cumbersome 
and it is useful to derive more transparent results in sim- 
ple limiting cases. At low temperatures the summation 
with respect to the Matsubara frequencies can be re- 
placed by the integration 27tT^ lj>0 — > J °° du>. In this 
limit we can obtain the analytical result for the average 
correction to the order parameter amplitude, Af^, 



A R 
ttT c 



d„^ 



U(A s0 /\A a0 \) 



Re[A a0 ] - A s0 
7bq|A q0 | 



(A12) 
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with 



U(a)= dz 



(z 2 + l) 3 / 2 Va 2 z 2 + l 



K(l-a 2 )- E(l-a 2 
T^a 2 



where K(m) — J*' 2 (1 — m sin 2 6) x l 2 dQ and E(m) = 
J* (1— msin 2 6) l t 2 d6 are the complete elliptic integrals. 
As we mentioned before, the uniform part of A s can be 
selected to be real. A,,.o = A^ . 

We can derive the simple analytical results for impor- 
tant particular case of (i) thin s-layer, d s -C £ s ,a> (h) 
weaker s-superconductor, A s o <C |A Q o|, and (hi) low 
temperatures, T <C T*. Due to the first condition, the 
dominating contribution to the gap correction is given 
by the coordinate independent part A S) 0) which is de- 
termined by the general formula (|A12I) . In the limit of 
A s q -C |A„o| we can use the asymptotics of the function 
U(a) in the limit a< 1, U(a) s» In {A/a) — 1, leading to 
the following simple result 



TlT r _ 



C y Re[A a0 ]-A s0 

d a 7Sa|A a0 | 



In 



4IA 



oO 



VsO 



-1 



(A13) 



The sign of A^ determines net effect of the s± supercon- 
ductor on the s superconductor, i.e., the sign of the prox- 
imity effect (positive vs negative proximity). The prox- 
imity is always negative in the TRSB state. In the case 
of the aligned state corresponding to Re[A a0 ] = A^o, 
as one can expect, s± gaps aligned with A s o enhance s- 
wave superconductivity while anti-aligned s± gaps sup- 
press superconductivity in the s superconductor. The rel- 
ative contributions are mostly determined by the electri- 
cal coupling between s-superconductor and the s± bands 
described by parameters 7g* ■ Another important factor 
is that the aligned gaps give positive contribution pro- 
portional to gap difference A a o ~ A s o, while antialigned 
gaps give negative contribution proportional to gap sum 
|A a o| + A s o- This give possibility to total negative prox- 
imity effect in the aligned state. 

Weak spatial dependence of A s (x) is determined by 
the components A S)OT with m > 0. At T = these 
components can be presented as 



A 



1 2 £U f\^*o\ \Re[A a0 ]-A s0 



7a A e* / ^ 



A s0 



iBa 



1 24 2 . 



zt m ^ 



(A14a) 

A^.U^O^el (A i4b) 



VsO 



iBc 



with (3 m = (■Km£ St A/d s y 

7 a • dz 

s,m 



{z 2 + 1)3/ 2 
1 



Vz^TT + p, 
| + V/3 2 l -im( v //^-i + /3 m ) 

l 



dz- 



(z 2 + 1) Vz 2 + S 2 V?+T + j3 

dz j3 m 



Vz^+l Vz^+l + Pm 
V Pra - 1 V J 

dz 



o Vz 2 +S 2 (Vi^TT + z?) 



In the limit /3 m >• 1 corresponding to d s <C £ S! asymp- 
totics of both Z a s m and Zf <m are zjffi « ln(2/3 m ). In 
the limits |A a0 | ^ A s0 and A s0 /3 m » |A q0 | correspond- 
ing to (3 >• 5 » 1, J a (5,f3) and J${8,fi) also have the 
same asymptotics 

Jo (z + /8) Vz 2 + <5 2 p o 
Collecting all terms, we obtain 

ln(A s0 /|A Q0 |) 



i+- 



ln 



2 {irm^s^/dsY 



2d s /C Aao-A s o 
(7rm) 2 7-Bq 



Using the relation (\x\— l) 2 = i+4 cos ( 7r " t:r ) we can 

approximately present the gap correction in real space as 

'(x+d s ) 2 1" 



A s (:r)«A s ,o--^— r- 

ln(A s o/|A a0 |) 



2d 2 



6 



2 In (7rV^.,A/da 



(A15) 



Correspondingly, for the Green's function in the same 
limits we derive 



$s(u),x) fa A s ,o +7rT t 



i + 



1 { x + df. 



E 



i 



2 V &,<. 



7sc 



(A16) 



In summary, simple analytical results given by Eqs. 
(|A13|) . (|Al5j) . and (|A16|) determine corrections to the 
s-wave gap and Green's function for a thin s-wave layer. 



2. s±-wave gaps and Green's functions 

We can evaluate corrections to the s± gap parameters 
and Green's functions following the same general route. 
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The difference is that the matrix structure of the self- 
consistency condition, Eq. (|5b|) . has to be properly ac- 
counted for. The first-order correction to $ Q with re- 
spect to the coupling strength 7^ is determined by the 
following equation and boundary conditions, 



f 2 <$>"_$ = _A 



G_s_ 

iBc 



■(A s0 - A q0 ), at x = 



(A17a) 
(A17b) 



and <!>' a = at x = d± with G 
D a /{2^ 



u/y/u* + \* a0 \ 2 ,e a ,u = 

|A q0 | 2 , and 



|A Q0 | 2 ) = e A |Aao|/V^ 
£,a A = D a /(2\A a o\). The self-consistency condition for 
corrections can be written as 

1 



UJ>0 



1\ 



A Q0 Re 




UJ 2 + 


A a0 | 2 



Aa 



^WapAp - In— A c 



(A18) 



with w a /3 = \~p — X~ 1 5 aj 3 and A is the largest eigen- 
value of the matrix X a /3- The matrix w a p is degenerate, 
WWW22 — W12W21 = 0, and its components are given by 



\W22 



A 2 _/4 + A 12 A 21 tA-/2 
detA ' 



W\2 =■ 



Al2 



detA 
(A19) 

with A_ = An ^ A22 and det A = A11A22 — A12A21. 

Similar to the s-wave case, we can split <fr a into the 
contributions induced by the boundary condition and by 
the correction to the gap parameter, <I> Q = $ Q ,b + & a ,A- 
The equation and the boundary condition for <f> at i,(x) are 



£ 2 $ 



n 

a.cj ^ a.b 



0. 



1 ^/u 2 + \A a0 \< 



The solution for § a ,b{x) is given by 



(A20a) 

(A s0 -A q0 ). (A20b) 



£, a>u y/uj 2 + \A a o\ 2 



S,a J Ba \Ju 2 + A 2 Q 

cosh [(a; — d±)/$ >a 



sinh(d±/^ Q 



(A s0 - A q0 ) 



(A21) 



The component <& a ,A(x) has to be found from the follow- 
ing equation and boundary conditions 

£,Aa " *«,A = -A a , (A22a) 
K,a = for x =0,d±. (A22b) 

We can again find <& a ,A( x ) and A Q (a;) using Fourier 
expansions, $ a< A(x) = J2 m $«,A,m cos (q m x), A a (x) = 



J2 m cos (q m x) with g„ 

we immediately find 



--rmt/d ± . FromEq. (|A22a.p 

rv . m 



- a. Zi.m 



> (A23) 

1 + e 2 ^g?„ 

As follows from the structure of the the self-consistency 
equation (| A18|) . the responses of the order parameters 
A a are different in the amplitude and phase channels. 
To proceed, we split all quantities into the amplitude and 
phase components, as illustrated in Fig. [2j X = X a +X^, 

xa = A a0 Re[XA* ] , ^ = x _ A a0 Re[AA* ] 



|a q0 



|a q0 



where X stands for A Q , <& a ,b, $ a ,A, and A s q. Such 
decomposition of A s0 is illustrated in Fig. [2] Explic- 
itly, we can write, A" = (Ai / 1 Ai 1) | A s0 1 cos</> and 
A s0 = -i(Aio/|Ai |)|A s0 | si n (ft. S ubstituting $ ajA into 
the self-consistency equation (IA18|) , we obtain the follow- 
ing equations for A° m and A^ m 



£ (wc* ~ K, m S a p) A%, m = 2nT £ 



ui>Q 



' (c 2 + |A q0 | 2 ) 3/2 ' 
(A24a 

1 



_(u J 2 + \A a0 \ 2 f /2 (l+e a , u ql) " 
*** 

Oi.b.m 



T ' 



[3 w>0 V w2 + l A a0h 



(A24b) 



uj>0 



v^ 2 +|A q0 | 2 i+£U<4 w 



-In 



where the Fourier components $^ 6 m and 6 m can be 
computed explicitly from Eq. (|A21I ) , 

K,b,m\ _ (2 - S m ) (d±&) 



7-Ba 



1 + £2. u « 



2 



y/^ + IA^p /Agp-A^ 
V^T+Al V A' 



sO 



s-0 



Solutions of Eqs. (|Al4l arc 



LJ 2( I> a 

AO _ o_T^ \ ^ j-ra /3,b,m 



At <m = 2^T ^ < Q/ 



( W 2 + |A /30 | 2 ) 3 / 2 ' 

/3,6,m 



where the matrices a/3 = \w a p 
two-band case are given by 



^ 2 + l^0| 2 



u m,al3 



W 22-^2 t ra 
-W21 



-W12 

Ea,< 
it 



(A25) 



(A26a) 



(A26b) 



the 



(A27) 
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The equation for the phase component at m — re- 
quires special attention. We note that the equation for 
the bulk gaps have the form 



E 



A^o = 



meaning that Eq. (|A24b|) at r. 
erate, i.e., its determinant vanishes, Dfff, — — rauE^ 



is actually degen- 



t»2Sf o + £? S 2,0 



'U,0 ~ -U7ii ^2,0 

Q. This degeneracy reflects gauge 
invariance with respect to identical phase change of the 
order parameters A a o. This means that the equation 
for A^ Q only has solution if its right-hand side satis- 
fies certain condition which, using the bulk gap ratio 

Aio/A 20 = -W 12 /(Wu - Ef )0 ) = -(w 2 2 - I^ i0 )/W21, 

can be written as 



2ttT 



W 2 lA W <Pf 



b.Q 



wi 2 A 2 o^ 



2,6,0 



uj>0 



0. (A28) 



The meaning of this condition is that the total "torque" 
from the interface forcing unform phase rotation of the 
s± gap parameters has to vanish. Using the relation 
W21/W12 — A21/A12 = vi/v 2 and Eq. (|A25|) we can 
rewrite this condition as 



Eq. (|A25[) for A^ determines interface-induced phase 
shifts (p a <C 1 of the averaged order parameters with re- 
spect to zero-order phases <j> and 4>—tt, see Fig. [2] These 
phase shifts are defined by relation A^ = iip a A a Q. The 
same phase shifts appear in the phenomenological frus- 
trated Josephson junction model. As the average phase 
shift (ipi + </?2)/2 can be absorbed into 0, we can set it 
zero and take tp 12 — ±<p/2. From Eq. (|A11|) . taking into 
account the above condition, we derive 



2ttT 



9 Z 



E 



I A, 



d±wi 2 |A 20 | ^ ^/uj 2 + A 2 sQ Ibi 



(A29) 



One can verify that this result does not change with the 
switching of indices 1 f> 2 in the right hand side. Using 
the expression for the partial Josephson energy, Eq. (|lll) , 
this result can be rewritten as 



E 



ip - 



J.l 



d±i>itwi2|A2o||Aio| 



(A30) 



sf) 



2t:t y 

with A^q = A s o sin <f>. Moreover, for the combination 
£a u v a/£, a we obtain 



D a V a 



1 



£a £ a y/uJ 2 + \A aQ \ 2 Ja^UJ 2 + |A q0 | 2 ' 

which allows us to present the condition in the form 



2 - T E ^ 



A q0 sin ( 



0. 



We immediately recognize that this condition is equiva- 
lent to Eq. (|A11I) (vanishing of the total Josephson cur- 
rent flowing through the interface). With this condition, 



For weak interband coupling, |Ai 2 |, | A21 1 <C An, A22, the 
parameter £12 = fiu>i2|A2o||Aio| represents the inter- 
band energy and this result coincides with the result ob- 
tained within the frustrated Josephson junction model 
in the case Eja -C d±£\ 2 - In this situation w\ 2 > and 
ip < 0. However, in contrast to this model, in a gen- 
eral situation the phases <p a do not fully determine the 
energy of the s± superconductor. Moreover, it was ar- 
gued that for the iron-based superconductors the pairing 
is dominated by the interband coupling, i.e., the oppo- 
site inequality holds, | A12 1 , IA21I An, X 2 2- In this case, 
as A12 < and A11A22 — Ai 2 A2i < 0, we have w\ 2 < 
meaning that ip > 0. Fig. factually illustrates this situ- 
ation. Using the results presented in this Appendix, we 
derive in the main text the amplitude of the sin 2(j> term 
in the Josephson current which determines the width of 
the TRSB state. 
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